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Abstract 

We discuss the hadroproductions of Ac, Ac,-D and D in the framework of the constituent quark- 
diquark cascade model taking into account the valence quark annihilation. The spectra of Ac and Ac 
in pA, TT A and ■n~ A collisions are well explained by the model using the values of parameters used 
in hadroproductions of D and D. It is shown that the role of valence diquark in the incident baryon 
is important for D productions as well as for Ac production. 

1 Introduction 

Hadroproductions of charmed particles have been measured in fixed-target experiments 012 and indicate 
the large leading/non-lcading asymmetries, i.e. asymmetries between production cross sections of leading 
particles and those of non-leading particles defined as 

^ fj(leading) — a{non — leading) 
cr(leading) + a{non — leading) 

The leading particle shares a valence quark in common with the incident hadron, while the non-leading 
one docs not. The leading particles are copiously produced at large Feynman x as compared with the 
non-leading particles in the forward region of the incident hadron. This is called as leading particle 
effect. For example, the asymmetry between D~{dc) and D^{cd) in ■K^{du) interaction with a nucleon is 
defined as A^- ^{D^ , D^) — {da{D^) ~ da{D^))/ {da{D^) + da{D^)) and increases from zero to nearly 
one with x in the tt^ fragmentation region. In the conventional pcrturbative QCD at leading order (LO), 
the factorization theorem predicts that c and c quarks are produced symmetrically and then fragment 
into D and D independently. As a consequence, LO pcrturbative QCD predicts no asymmetry in contrast 
with the experimental data. The next to leading order (NLO) calculations generates asymmetries but 
they are much smaller than the data^ ^ . 

To explain the leading/non- leading asymmetry, many attempts have been investigated: fcr factorization^, 
string fragmentation0f7','S^, intrinsic charm contributions'^, recombination process with beam remnants |10l 
[m E. 12* Ji» US] , recombinat ion with participants in the hard scattering process [THl 1 17j . recombination 
with surrounding light quarks in projectile and target [THj. recombination using valon concept [T^ . light 
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quark fragmentationpilj . meson cloud model^T) and so on. Productions of Ac and Ac in T,~A |21I221 and 
T,~A,pA and tt~A collisions|23 have been measured and analyzed in Refs.^lElElEl and j^S]. These 
models explain the leading/non-leading asymmetry of charmed hadrons successfully, however, there are 
some uncertainties in the perturbative QCD calculations as reviewed in Refs.|26) and |27j . Hadronization 
of heavy quarks is still an open problem, since it is intrinsically non-perturbative process. 

We have applied the covariant quark-diquark cascade model for hadroproductions of charmed mesons 
as well as hadrons composed u,d and s quarks [2HI UHl EOl E]- In the present paper, we investigate the 
hadroproductions of Ac, Ac,!) and D and the leading/non-leading asymmetry in the framework of the 
constituent quark-diquark cascade picture and discuss the role of constituent diquark. 



2 The model 

2.1 Valence quark distributions in the incident hadron 

We assume that a baryon is composed of a constituent quark and a constituent diquark, and a meson 
is composed of a constituent quark and a constituent anti-quark. The constituent quark (diquark) is 
a quark-gluon cluster consisting of the valence quark (diquark), sea quarks and gluons. Hereafter, we 
simply call them as quark and diquark. When the collision between the incident hadrons A and B occurs, 
it is assumed that the incident hadrons either break up into two constituents or emit gluons followed by a 
quark-antiquark pair creation. There are four interaction types: a) non-diffractive dissociation, b) and c) 
single-diffractive dissociations of B and A, and d) double-diffractive dissociation types as shown in FigQ] 
The probabilities of these types to occur are (1 — Pgi)'^,Pgi{l — Pgi),Pgi{l — Pgi) and P^;, respectively. 
Here we denote the quark-antiquark pair emitted from A [B) via gluons as Ma ( Mb)- The probabilities 
of Ma ( Mb) to be uu, dd, ss and cc are are denoted by Puu-, Pddi Pss and Pee, respectively. 




Figure 1: The interaction mechanism in AB collision: (a) Non-diffractive dissociation type, (b), (c) 
Single-diffractive and (d) double diffractive dissociation type mechanisms. 

The light-like variables of A and B are defined as follows: 

^ _ ± Pcm B _ T Pcm 

^a± — 1= , ^o± — 1= , 1^7 



where ^/so = E^ + E^ is the center of mass energy of the incident hadrons A and B. The momentum 
fractions of A,B and Ma for the process shown in Fig. lb are given as follows: The momentum fraction 
of Ma is fixed by the distribution function 

HM^/Aiz) = z^^'-\l - z)f''--'/BiPgi,pM), (3) 
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and the uniform distribution R in the interval from zero to one as, 

— Xq_^z^ x_ — Xq R, (4) 

Then the momentum fractions of the incident particles A and B become as follows: 

x^^x^+{l-z), x^=ml/{xi so), 

x_ = Xq_ — {x_ — Xq_{1 — R)), x^ = Xq^, (5) 

where the mass shell condition is considered and transverse momenta are neglected. Analogous expressions 
describe the momentum fractions of A, B and Mb for the process shown in Fig.lc , upon replacing 
superscripts A B and subscripts + <-!■—. The single-diffractive dissociation type mechanisms result in 
diffractive peaks of spectra of the same kind of the incident particles A and -B at |a;| « 1. The probability 
of the incident hadron to emit gluons followed by a quark-antiquark pair Pgi and the parameters figi and 
Pid in ® are chosen to reproduce the diffractive peaks of the spectra such as ■k'^,K'^ and p in ■k'^p, K'^p 
and pp collisions at \x\ « 1^^. The parameters are chosen as Pgi = 0.15, /3g; = 0.1 and — 3.0. 
However, the main contributions to charmed hadrons come from the non-diffractive dissociation type 
mechanism. 

The distribution functions of the constituents in the projectile A composed of a and a' are described 

as 



z 



Pa 



Similarly, the distribution function oi q in Ma is given as iJ^/M^l-z) = z^''~^{l — z)^''~^ / B{[3q, (iq). The 
dynamical parameters /3's in JB)), which determine the momentum sharing of the constituents, are related 
to the intercepts of the Regge trajectories as /3u = /3(j = 1 — o-p-uj{^) ~ 0.5, /3s = 1 — a^{Q) ~ 1.0, /3c = 
1 — aj/^(0) ~ 4.0 |32l For flavor anti-symmetric diquark [qq'] and symmetric diquark {<z<z'}, the 
dynamical parameters (3[qqi] and P{qqi] are chosen as I3[qqi] = 7[ ](/3g + Pq') and /3{gq'} = 7{}(/3g + Pq'), 
respectively. The commonly shared constituents of Ac with the incident p, and ■n~ beams are ud- 
diquark, d and d-quark respectively. Ac has a commonly shared u quark with ■n~ and K~ beams and Ac is 
a non-leading particle in K~ beam. These distributions of the constituents in the incident hadrons affect 
the shapes of spectra of produced hadrons significantly. Hadrons are produced on the chain between 
a constituent valence quark((j'^)(anti-quark(g^)) from the incident hadron A{Ma) and the constituent 
valence diquark((7(7B(gs))(quark((7B)) from the incident hadron B{Mb) for a meson-baryon collision as 
illustrated in Fig.l. 

The small intrinsic transverse momenta of constituents in the incident hadrons are introduced by 
changing the direction of the incident constituents by an angle in the rest frame of the chain. Here 
we choose the distributions for z = cos ipin as 

DM = ^0^{l + -A'- (7) 

in the region — 1 < z < 1. Here we choose very small intrinsic transverse momenta and use the value 
Pin = 40 used in Ref.[3J in order to reproduce p^ dependences of D mesons. 

We take account of the annihilation of the valence antiquark from beam and valence quark from target 
in a slightly different way from the one pointed out in Refs. ^Hl and JI]. When the incident constituents 
in A{Ma) and B{Mb) are quark and its anti-quark, it is assumed that the annihilation process occurs 
with the probability Pan and the non-annihilation process occurs with 1 — Pan- Let us consider the case 
of TT^p collision. The annihilation process 

u + u—>q + q^ "hadrons" + "hadrons" (8) 
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occurs with the probabihty Pan, when u in 7r~ interacts with u in p. Branching ratios of uu — > uu, dd, ss, 
and cc are chosen to be equal to Pun, Pdd, Pss and Pee for the channels allowed energetically. The produced 
q and g in © are supposed to be non-free in terms of hadronization mechanism. Thus it is assumed that 
q has tendency to be produced in the forward direction of u in 7t~ beam. Therefore, here we choose the 
distributions for z — cos (pan as 

DU^) = ^-tlil + zf^^ (9) 

in the region — 1 < z < 1, where fan is the angle between the directions of u and q. To reproduce 
the yl^-jv(-D*', -D"), we choose the value Pan — 0.2 and use a value of Pan — 20 which roughly means 

Daniz) ^ S{Z - 1) 31,. 



2.2 Quark-diquark cascade process 

Hadrons arc produced on the cascade chain by processes 



q ^ Miqq')+q', (10) 

q_ ^ B{q[q'q"]) + [q'q"], B{q{q'q"}) + {q'q"}., (11) 

[q'q"] ^ B{q[q'q"])+q, (12) 

[q'q"] ^ M{ qq') + [qq%M{qq') + {qq"}, (13) 

{q'q"} ^ B{q{q'q"})^, (14) 

{q'q"} ^ M{qq') + [qq%M{qq') + {qq"}, (15) 



where q denotes u, d, s and c and [q'q"] does [ud], [us], [uc], [ds], [dc], [sc] and so on. The symbols [ ] and { } 
denote the flavor anti-symmetric and symmetric diquarks, respectively. Meson production probabilities 
from q, [q'q"] and {q'q"} are 1 — e, 77[ j and rj^ j, respectively. The parameter e is the probability of baryon 
productions from a quark and is related to the baryon productions from tt and K beams resulting in 
e = 0.07. The processes (|13|) and (I15II contribute to production of baryons which share no valence 
quark or diquark in common with the incident hadron. Although both and Q share no valence quark 
common with the incident proton, the spectrum of is harder than that of 51 in the proton fragmentation 
region. In our model this is described as follows: after emitting K mesons through these processes, the 
incident diquark in the incident proton converts into {ss} diquark and produce Q by the process H14|l . 
The difference between the spectra of 57 and Ct in pp collision is described by the processes (|13|) . (|15ll and 
(I14|l and the parameters are chosen as ?7[] = r]^y — 0.25. |28l IHlj 

For the cascade process a H{ab) -I- & in the beam (target) side, the distribution function of light-like 
fraction of b is assumed as 

Ftaiz) = Z^^^-\l - zf^+P^-^lB{-lPa,Pa+Pb) (16) 

and the light-like fraction of h becomes = zx'i^ (x^L = zx"i_). We fix the transverse momentum of the 
hadron _ff as a difference between two random points = p-^j ~ Pti with their lengths fixed by the 
distribution function 

G(p|) = ^exp(--^p|) (17) 

in space, where niH denotes the mass of H . The parameter C is fixed from the experimental data 
on distributions of pionspU]. We use the values of 7 = 1.75 and C = 1.8GeV~5, from previous 
analyses [3D| 12] ■ 
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We now see how the cascade processes occurs between the incident constituents, for example, q from 
A (Ma) and {q'q"} from B. We redefine the hght-hke fractions of q and {q'q"} in the rest frame of 
the incident constituents q and {q'q"}. The momentum sharing of the cascade process q + {q'q"} 

H{qb) + b + {q'q"} from q with and {q'q"} with x^^'/'^ takes place as follows |21|33|: Using Ea. ((TC|l 
with a = q, wc fix the light-like fractions of b and H{qb) as = x'^^z and x^ = — x'!^_, respectively 
and put x'L — x'L- The transverse momentum of H is fixed by Eq. Ijl?!) . Then x^ is fixed from the 
onshell condition x^x^ — {rnjj +Px)/s', where ^/s' is the subenergy of the incident q and {q'q"} system. 
The transverse momentum of b is = p\. — p^. The light-like fraction of {q'q"} is decreased to 
j{9 9 } _ 9 } _ r^ii_ jf energy of the b -f {q'g"} system is enough to create another hadron, the 
cascade process such as 6 -|- {q'g"} ^ & + -ff (c{(7'g"}) -I- c takes place in the opposite side. In the final 
step, we assume that the constituents recombine into one or two hadrons according to the processes: 
q + q' ^ M{qq'), q+[q'q"] ^ B{q[q'q"]), {qq'} + [q"q"'] M {(^) + M (q'q^) and so on. The momenta 
of the recombined hadrons are the sum of those of the final constituents and are offshell. In the cms of 
incident hadrons A and B, we have energy- momentum conservation relations = and — ^/sq, 

i i 

where Pj denotes the three momentum of «-th produced particle. In order to put recombined hadrons 
onshell, we multiply the three momenta of all produced hadrons by a common factor / so that the 
summation ^ \/ (/pj^ + "^f would be equal to 

i 

We take into account pseudoscalar (Ps), vector iV) and tensor (T) mesons. The production probabil- 
ities for them are assumed to be Pps = 1/9, Py = 3/9 and Pt = 5/9. For baryons we consider lower lying 
baryons: octet (iJg) and decuplet (Bio) baryons composed of u,d, and s fiavors and the corresponding 
ones with charm flavor. Octet and decuplet baryons are described as 

IBs >= COse\q[q'q"] > + sin 9\q{q' q" } >, (18) 

\Bio >= \q{q'q"} > ■ (19) 

We assume the production probabilities for octet and decuplet baryons from symmetric diquarks to be 
Pbs = 1/3 and Pbw = 1 — Pbs = 2/3, respectively. Then the emission probabilities of individual cas- 
cade processes are determined from the above probabilities. For examples, the probabilities to produce 
7r+, p+, aj,..., A++, p, A+, S++ and S*++ from a u quark are as follows: 

(1 - e)P,3Pps, (1 - e)P^,jPv, (1 - e) P,jPt,..., 
e Pj, e P™ P^rf- i Pb8 sin^ 9/{ i Pbs sin^ + | Pgio) , 
e Pun Pdd I Pbw/ {\ Pbs sin^ + | Psio),-, 
e Pj f Pbs sin^e/ (| Pbs sin^ + i Pbw ), 
e Pj i PbwI (1 Pbs sm^ + \ Pbw), 

where the factors 1/3 and 1/2 are flavor SU(4) factors. Directly produced resonances decay into stable 
particles. 

The mixing angle of symmetric and anti-symmetric diquarks in octet baryon cos and the pair creation 
probabilities of uu, dd and ss are determined from the analyses of non-charmed hadron productions in 
irp, Kp and pp collisions. In Ref.^Hli we used cos^ 9 = | and Puu = 0.45, P^d — 0.45 and Pss = 0.10. 
To reproduce the data on ir^p D^X cross section, we choose cos^ 6* — \,Pcc = 0.0003, Puu — 
0.435, Pdd = 0.435 and Pss = 0.1297|2I]. The dynamical parameters 7[ j and 7{} as well as the mixing 
angle of symmetric and anti-symmetric diquarks in octet baryon is related to the tt^ spectra in proton 
fragmentation region and are chosen as 7[ j = 1.5 and 7{j = 2.0 |3U[I5T] . 
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Figure 2: Fcynman a; distributions of A+ and in a) p, b) E , c) tt , d) iiT and charmed mesons 
distributions in e) and f) collisions with N and Pb targets. 



2.3 Interaction with nucleus 

We now consider hadron-nucleus interaction. We regard the hadron-nucleus collision as a sum of hadron- 
nucleon collisions and neglect the intra-nuclear cascade of produced hadrons. The projectile hadron 
successively collides with nucleons inside the nucleus. The probability PhAii^) that the incident hadron 
h collides with nucleons is calculated by the Glauber- type multiple collision model[2n|- The number v 
is determined by the distribution of nucleons in the nucleus and the cross section of the incident hadron 
with a nucleon ahN- We take the values apN = 40, cr^-Af = 34,cr^-jv = 24 and (Tx-n = 21mb in the 
hadron-nucleus collisions. The nucleon number density of the nucleus with a mass number A is assumed 
as 

PM =POT— ^^73f7T, (20) 

l + exp(— j^) 

where we choose ta — l.l9Ai — l.QlA~i (fm) and d = 0.54 (fm). Then it is assumed that AdiNi, M2N2, 
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Mi,-iN,y-i and hNi, interactions take place, where AU is qq state emitted from h as Ma discussed in 
section At each coUision, the projectile hadron h looses its momentum. From the data on A depen- 
dences of the spectra of the projectile hadrons, we put the degradation rate of the momentum of h as 
P(z) = 0.25z"-25-i which is used in Ref.,29,. 

3 Comparison with experiment 

In this section, we give the results of our model for charmed hadron productions in p, E^, tt^ and 
beam interactions with nuclei. To see the target mass number (A) dependence, we show the x distributions 
of Ac,Ac,D and D productions in p, I]~,7r~ and K~ interactions with N and Pb targets at = 600 
GeV/c in Fig. 2. All spectra increase with A in the region x ^ Q. The spectra of become soft with 
A in the baryon beam fragmentation region. This is due to the momentum degradation of the incident 
baryon through the multiple interactions with nucleons in the target. 

The spectra of Ac scarcely depend on A in the baryon beam fragmentation regions. In meson beam 
fragmentation regions, both Ac and Ac slightly depend on A and become soft with increasing A. In 
baryon beam fragmentation regions, A-dependences of the spectra of D mesons are smaller than those 
for D mesons. The shapes of the spectra for D mesons in baryon beam fragmentation regions are steep as 
compared with those of D mesons in baryon and D and D mesons in meson beam fragmentation regions. 

In Fig. 3, we show the x distributions of Ac and Ac in p, S", 7r~ and K~ beam interactions with C and 
Cu targets in arbitrary unit and compare the results with the experimental data E78l[^. The shapes 
of Ac and Ac in p, and fragmentation regions are well explained by our model. In our model, 
the shapes of the spectra are strongly related to the distributions of valence constituents in the incident 
hadrons, which are characterized by the parameters /3's in 10. Diquarks in the incident baryons have a 
tendency to bring a large momentum fraction of the incident baryons. The valence constituents d and u 
in n~ beam have the same distribution i.e. energetic d and wee u or vice versa. The valence constituents 
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Figure 4: Asymmetries between A+ and in a) p b) S , c) tt and d) collisions with C and Cu 
targets. The experimental data are taken from Ref. 21. 



d and u in 7r~ are harder than the valence quarks in and p beams. The ds diquark in is harder 
than the ud diquark in p. The valence {ud} diquark in p beam is the hardest among the commonly 
shared constituents between Ac and the incident p, S~,7r^ and K~ beams. Therefore x-distribution of 
Ac in proton beam is harder than those in other hadron beams. 

The asymmetries between the spectra Ac and Ac in pA, Y^~A, 7t~A and K~ A collisions are shown in 
Fig. 4. The results are in good agreement with the experimental data except for 7r~ beam. For 7r~ beam, 
Ac and Ac share the common valence constituents d and m, respectively. Therefore the same spectra are 
expected for Ac and Ac productions in the region of 7r~ beam fragmentation. In our model, however, the 
excess of Ac over Ac is seen at 0.3 ^ x. This is due to the fact that energetic Ac and Ac are produced in 
the c?-quark-diquark and the fZ-antiquark-quark chains, respectively. The former chain is slightly shifted 
to the target fragmentation region as compared with the latter one due to the large momentum fraction of 
the diquarkjSD]. Therefore our model gives the negative asymmetry Att-aCAc, Ac) at 0.3 ^ x. The valence 
diquark in S~ beam contributes to Ac production via a processes such as ds — > + K~ c + Ac -I- K~ 
in fragmentation region. The contributions to Ac production from the valence diquarks dominate the 
contributions from the valence quarks in E~ beam. 

In Fig. 5, We show the results of asymmetries between charmed and anti-charmed mesons in p and 
interactions with N,C and Cu targets. The asymmetries ApA{D^ , D^), ApA{D^ , D'~'), A^^- ^{D^ , Df) 
and A-^-a{D^,D^) are small as compared with the asymmetries Ap^(Ac,Ac) and Ay^-a{Ac,Ac). Al- 
though both Z?+ and D~ in p and and in E~ beam are non-leading particles, the results of 
asymmetries ApN{D^ , Df) and A-^- ^{D'-' , D'^) are non-zero and have large values. This is because in 
our model D mesons are produced mainly from the chain between the valence diquark in the incident 
baryon and the valence quark in the target nucleon. For example, Dj in proton fragmentation region 
is produced via a process such as ud sd + cA + + . The more energetic the valence 

diquark becomes, the more copiously D mesons are produced in the baryon fragmentation region. In our 
model the valence ds diquark in beam has a larger momentum fraction than the valence diquarks in 
p beam. Consequently the asymmetries A-^- ^{D'^ , D'^) are large as compared with ApA{D'^ , Df). The 
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Figure 5: Asymmetries A{D, D) between leading and non-leading charmed mesons in a) pA and b) Yi~A 
collisions. Asymmetries between non-leading charmed mesons in c) pA and d) collisions. The 

experimental data are taken from Ref. 20. 

asymmetry A-^- j^{Dj , Df) is larger than other charmed meson asymmetries in p and E~ beams. This 
is also due to the large momentum fraction of the valence ds diquark in E~. In the target fragmentation 
region, both D and D mesons increase similarly with the mass number A and the asymmetries show small 
mass number dependences. Since the spectra of D mesons increase with A and those of D mesons change 
little, there are considerable target mass number dependence in the asymmetries A{D, D) at ^ a;. 

4 Conclusions 

1. We have examined hadroproductions of charmed baryons in the frame work of the quark-diquark 
cascade model. Their spectra are well explained by using the values of the dynamical parameters fixed 
from Regge intercepts. The large leading/non-leading asymmetries ApA{-^c:^c) and ylg-^(Ac,Ac) are 
naturally explained. 

2. The valence diquark in the incident baryon, which brings a large momentum fraction of the beam, 
plays an important role for D productions as well as for Ac Droduction|28j. The non- leading anti-charmed 
meson productions in baryon beams are largely affected by the valence diquark distribution in the incident 
baryons through the processes (13) and/or (15). 

3. In meson beam fragmentation regions, A-dependences of the spectra of non-leading charmed mesons 
are small as compared with those for leading charmed mesons. There are considerable A-dependences in 
Ac and D meson productions, while those for Ac and D productions are small in baryon beam fragmen- 
tation regions. 

4. There are discrepancies between the data and calculations of asymmetries Aj^- a{Kc,Kc) at 7r~ 
fragmentation region as seen above. Piskounova discussed this problem successfully by introducing the 
transfer of string junction in quark-gluon string model|3S]. To explain the experimental data in our 
model, it is necessary, for example, a energetic d quark from 7r~ beam and a wee uc diquark from target 
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nucleoli to recombine into a A^. Here the uc diquark is converted from the valence diquark in the target 
nucleon by emitting several mesons through the processes (13) and/or (15). The significant value of the 
asymmetry A^-j^{Ac,Ac) as well as the small value of ^4^-^(13", Z)") in meson beam fragmentation is a 
challenging problem. 
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